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Abstract: The geomagnetic field causes not only the East-West effect on the primary cosmic rays but also affects the 
trajectories of the secondary charged particles in the shower, causing their lateral distribution to be stretched along certain 
directions. Thus both the density of the secondaries near the shower axis and the trigger efficiency of a detector array 
decrease. The effect depends on the age and on the direction of the showers, thus involving the measured azimuthal distri- 
bution. Here the non-uniformity of the azimuthal distribution of the reconstructed events in the ARGO-YBJ experiment 
is deeply investigated for different zenith angles on the light of this effect. The influence of the geomagnetic field as well 
as geometric effects are studied by means of a Monte Carlo simulation. 
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Cosmic rays (CRs) are charged particles and their paths are 
deflected by the magnetic fields. The galactic magnetic 
field randomizes the CR arrival directions. The geomag- 
netic field restrains low-rigidity CRs from reaching the ter- 
restrial atmosphere and causes that the CR flux is lower 
from East than from West. The geomagnetic field acts also 
on the charged particles of the extensive air showers (EAS) 
during their path of few kilometers in the atmosphere. Coc- 
coni [ 1 1 suggested that the lateral displacement induced by 
the Earth magnetic field is not negligible with respect to 
the Coulomb scattering when the shower is young. There- 
fore the effect could increase for high altitude observations. 
In particular the shower extension along the East- West di- 
rection is larger than along the North-South direction (the 
opposite for the particle density). The different density of 
charged particles introduces an azimuthal modulation due 
to the different trigger efficiency of EAS detectors. This 
modulation was observed at the Yakutsk array for EAS with 
energy above 50 PeV [2 1 and at the Alborz observatory for 
energy above 100 TeV J3). A North-South asymmetry has 
been observed also in EAS radio-experiments [4 1 and it can 
be explained as an effect of the geomagnetic Lorentz force 
on the EAS charged particles. Therefore the geomagnetic 
influence on the lateral distribution must be taken into ac- 
count in the EAS simulation [5 1 and it was also suggested 
that a pointing correction is necessary for Cerenkov tele- 
scopes because of the geomagnetic field |6|. 



The effect of the geomagnetic field on the trigger effi- 
ciency of the ARGO-YBJ experiment was studied and sim- 
ulated [71- Due to the field B the average shift (d) of an 
electron (or positron) in the shower plane is 



qh 2 Bsinx 
2E P cos 2 6 



(1) 



where \ is the angle between B and the particle velocity 
v, q the electric charge, h the average vertical height of the 
electron path, E e the average energy and the zenith angle. 
The trigger efficiency is connected to the size of the EAS 
footprint (lower particle density corresponds to lower trig- 
ger efficiency). Therefore the stretching of the EAS lateral 
distribution introduces a modulation with respect to the az- 
imuth angle ((f)). The azimuth distribution is expected to 
be 



dN 



= N {l+ gi c 



+ <t>i)+g 2 cos [2 (<£ + <fe)]} 



(2) 

with <pi = <j>2 — — 4>b (4>b is the azimuth of the geomag- 
netic field), gi oc sin26 and <?2 oc sin 2 6 if the modulation 
is totally geomagnetic. Furthermore Q 
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siwub sin 2 9 
where 9b is the zenith angle of the geomagnetic field. 
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Figure 1 : Absolute value of sin\ versus local angular co- 
ordinates (8 and (f>) of the arrival direction of a charged par- 
ticle in the ARGO-YBJ reference frame, x i s me angle 
between the magnetic field at YangBaJing and the velocity 
of the particle. 

1 ARGO-YBJ experiment 

ARGO-YBJ is a full-coverage EAS experiment, located 
close to the YangBaJing village in Tibet (People's Repub- 
lic of China) at 4300 m above sea level. Its geographical 
coordinates are 90°31'50"£: and 30°06'38"N. The exper- 
iment is mainly devoted to Very High Energy (VHE) 7- 
astronomy and CR studies. The detector is essentially a 
continuous carpet (78 x 74 m 2 ) of Resistive Plate Coun- 
ters. The detection area is enlarged to 110 x 100 m 2 by 
means of a partially equipped guard ring. The time-space 
pattern allows a detailed reconstruction of showers induced 
by gamma and charged primaries. 

In the ARGO-YBJ reference system the azimuth angle of 
EAS is defined with respect to the detector axes in the an- 
ticlockwise direction (<f> = 0° for showers aligned with the 
x-axis and moving towards the positive direction). Thus the 
azimuth angle of showers going towards the geographical 
North is 4> N = 71.96° ± 0.02°. The quoted error is due 
to the measurement of the orientation of the detector axes 
with respect to the geographical reference system. 

According to the International Geomagnetic Reference 
Field (IGRF) model available on the NOAA web site 
at YangBaJing the geomagnetic field (B = 49.7 fj,T) has 
the following angular coordinates in the ARGO-YBJ refer- 
ence system 

9 B = 46.4°, (j) B = 71.89°. (4) 

The geomagnetic effect on primary CRs is negligible for 
the EAS collected by ARGO-YBJ. Thus we focus on the 
effects on the secondary particles in the shower. The 
Lorentz force acting on charged particles depends on the 
angle \ used in Eq. (Q~|). The absolute value of sin\ is 
shown in Fig. Q] as a function of 9 and </>. It has been 
checked that the variation with the height is negligible. In 
fact from 4.3 km up to 30 km above sea level the product 
Bsin\ changes less than 2%. 



2 Data analysis 

This analysis is based on the data collected in 3 days (Octo- 
ber 12-14, 2010) with the trigger condition that at least 20 
time-pixels (pads) are fired. In order to get a reliable recon- 
struction of the shower direction, the following cuts have 
been applied: shower core reconstructed inside a square of 
40 x 40 m 2 at the center of the carpet, zenith angle lower 
than 60°. After these cuts ~ 130 millions of events are 
selected. 

The timing calibration of the pads has been performed ac- 
cording to the Characteristic Plane method (9)- This proce- 
dure removes the systematic time differences and arranges 
the time correction in order to make null the mean value 
of the direction cosines. After this step a last premodu- 
lation correction is necessary in order to remove a small 
over-correction intrinsic to the method [10|. It is remark- 
able that small systematics in the timing calibration prevent 
the azimuthal analysis and the premodulation correction is 
crucial in order to get the proper (^-distribution. 

The azimuthal distribution is shown in Fig.|2]together with 
a fit performed according to function @. x 2 l n df ~ 4 
and this is due to a small inefficiency at <f> ~ n X 90° (n = 
0,1,2, 3, 4) which has not been examined here but does not 
invalidate this analysis. The fit results are 

,9i = (1.521±0.012)%, 0i = -72.43°±0.47° ~ -0s, 

.92 = (0.587±0.012)%, 02 = -86.55°±0.60° ^ -<j) B - 

The first coefficient is 3 times higher than the second one. 
The phase of the first harmonic results fully compatible 
with what expected. This is not the case for the phase of 
the second harmonic. 

The geomagnetic origin of the modulation can be checked 
also by studying the dependence of the harmonic coeffi- 
cients on the zenith angle. According to Q 

.91 = k\ sin29, 92 = ki sin 2 9. (5) 

This behaviour is verified for gi (see the fit in Fig. [3]), not 
for 32 (the plot is missing for room saving). 
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Figure 2: Azimuthal distribution of EAS sample. Fit with 
function (0 is superimposed. 
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Figure 3: Coefficient g\ versus the zenith angle. Function 
91 = k\ sin29 is used for the superimposed fit. 

The unexpected result for the second harmonic cannot be 
explained as a consequence of some systematic error. The 
disagreement can be solved simply assuming a mix of geo- 
magnetic and detector effects. The periodicity of a detector 
effect must be 180° and then acts only on the second har- 
monic, its phase must be 0° or ±90° in connection with x 
and y-axes. The fit of Fig. [2] suggests the second solution 
because <\> 2 = —86.55° is in-between —4>b and —90°. 

In this hypothesis the second harmonic can be splitted in 
two parts: one (2B) is due to the magnetic field, the other 
one (2 A) can be originated by a detector asymmetry. Three 
different data sets have been selected on the basis of the 
zenith value in order to disentagle these two effects. The 
^-distributions of the 3 subsamples (a, (3 and 7 in Fig. 
can be fitted all together with a single function: 




= iV 4 {1 + ki{sin29)i cos (0 + ^1) 



+ k 2B {sin 2 e)i cos [2 O + 0i)] 

+ g i 2A cos[2(<f> + (f) 2A )}} (6) 

where the coefficients of the magnetic component are writ- 
ten using Eq.s ((5]), the magnetic phase is the same for 
first and second harmonic and the index i = a, f3, 7 in- 
dicates the subsamples. The quantities Nu (sin26)i and 
(sin 2 0)i have been estimated separately for each subsam- 
ple. Then the fit parameters are k\, k 2 s, 4>i> 9 2 a> 92A> 91a 
and 4>2A- The result of this new fit is reported in Tab. [TJ 
The value of k\ is compatible with the fit in Fig. [3] The 
ratio k\jk 2 B = 2.41 ± 0.90 is compatible with 1.9, value 
calculated according to (0 and (0]). The magnetic phase is 
in perfect agreement with the expectation {(pi = — </>_b). It 
is remarkable that g l 2A increases with 8 and c/) 2 a is compat- 
ible with —90° as expected for a detector effect. 

A possible explanation of the detector origin of the har- 
monic 2 A is that the pads have a different density along the 
detector axes (1.54 pads / m along x-axis and 1 . 76 pads /m 
along y-axis). Taking into account the trigger requirement 
(at least 20 pads fired in a time-window of 420 ns) the dif- 
ferent density can explain the higher trigger efficiency for 
showers along the y-axis (<f> = ±90°). The increasing of 



ki (%) 


2.101 ±0.016 


k 2B (%) 


0.87 ±0.32 


(°) 


-71.86 ±0.44 


9%A (%) 


0.142 ±0.024 


9 & 2A (%) 


0.350 ±0.061 


9l A (%) 


1.25 ±0.13 


4>2A (°) 


-92.5 ±2.9 


X 2 /ndf 


491/209 



Table 1 : Results of the fit with function (0 of the distribu- 
tions in Fig. [4] 



the effect with the zenith angle is the expected consequence 
of this hypothesis. 

In order to check the results for the first harmonic we ob- 
serve according to function (0 that the mean values of the 
direction cosines (I) and (m) of the showers depend on k\ 
and </>i as in the following 

k\ 

(I) = +—cos(f>i(sin29 sinO), 
k\ 

(m) = — —sin(f>i(sin29 sinO). 

In the analyzed sample the mean values are 

(l) = (12.12 ±0.29) x 10~ 4 

(m) = (35.56 ±0.30) x 10~ 4 
{sin 29 sin 9) = (3588.39 ± 0.20) x 10~ 4 
As a consequence 

h = (2.094 ± 0.016)%, 0i = -71.18° ± 0.57° 

These values are fully compatible with the the fit results in 
Tab.Q] 

3 Simulation 

Beams - In order to check the effect of the magnetic 
field on the data collection and on the EAS reconstruc- 
tion, some "beams" of primary protons have been simu- 
lated with 9 = 45°, E = 3 TeV and interacting at 19 km 
of height. The Corsika code ifTTI has been used to simulate 
the shower development and a GEANT-based code ifTSll to 
reproduce the detector response. Nine different beams have 
been generated with 3 different values of the azimuth an- 
gle ((f) = 71.5° where sin\ ~ 0, cj> = 251.5° where sin\ 
is maximum and the intermediate angle <f> = 161.5°) and 
assuming 3 different values of the magnetic field (null, the 
effective magnetic field and twice the effective magnetic 
field). Negative and positive EAS components have been 
studied separately. It has been verified that the geomag- 
netic field separates positive and negative cores of about 

4 m right on the East- West direction. This stretching of the 
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Figure 4: Azimuthal distribution in different 0-ranges (0° - 20°, 20° - 40° and 40° - 60°). Fit with function © is 
superimposed. 



(f> sin X B = fiT 49.7 fjT 99.4 fjT 

71.5° 0.02 29.3% 29.4% 29.3% 

161.5° 0.87 29.6% 28.4% 26.0% 

251.5° 1.00 29.0% 27.7% 24.9% 



Table 2: Trigger efficiency for simulated "beams" of CRs 
(see Sec. 0). Statistical error is 0.1%. 



lateral distribution does not affect the combined reconstruc- 
tion (positive + negative component) but acts on the trigger 
efficiency, according to [7 |. The results are summarized in 
Tab. |2] and confirm that the trigger efficiency decreases as 
the geomagnetic field or the sinx value increase. 

Complete sample - A sample of proton and helium- 
induced EAS has been simulated in the energy range 
10 GeV — 10 PeV, with a uniform azimuthal distribution 
and the proper geomagnetic field. The analysis chain of the 
real data has been applied to the simulated ones, with the 
exception of the timing calibration that is not necessary for 
simulated data. The azimuthal distribution of triggered and 
selected events (~ 2 x 10 6 ) is not uniform (Fig. [5}- The 
shape is similar to that of Fig. [2] and the parameters of the 
fit according to function (0 are almost compatible (the co- 
efficients show larger discrepancies) but with large errors 
of the simulation fit. A larger simulation will be necessary 
in order to fix these few percent effects. 

4 Conclusions 

The modulation of the azimuthal distribution of a large 
EAS sample has been analyzed. The origin and the fea- 
tures of this modulation are fully understood. It is well 
described by means of two harmonics, the first one of the 
order of 1.5%, the second one of the order of 0.5%. The 
first harmonic is due to the geomagnetic Lorentz force on 
the shower charged particles. The second harmonic is the 
sum of magnetic and detector effects. 

It is remarkable that in addition to the Moon shadow anal- 
ysis the absolute pointing accuracy of EAS arrays could 
be tested also with deep studies concerning the azimuthal 
modulation. 
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Figure 5: Azimuthal distribution for simulated data and fit 
with function (fJJ. The statistics does not allow to use func- 
tion ©. 



